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1. Introduction 

\Q ■ 

o : 

[-x. , The study of the Laplacian on a metric graph, a concept known as quantum graphs, 

now serves as a toy model for quantum chaos [H [2j [3]. Indeed, there exists an exact 

.£h ! trace formula relating eigenvalues and periodic orbits [U El H]. Moreover, depending 
on the graphs, exact computations of these orbits may be possible, whereas they are 
out of reach in most dynamical systems. It has thus been shown PQ El E] that spectral 
statistics of simple generic graphs follow random matrix statistics when the size of the 
graph tends to infinity, as expected for chaotic quantum systems. Star graphs (graphs 
formed by a central vertex connected to v other vertices by edges of different lengths) 
play a special role because of the high degeneracy of their periodic orbits. As could be 
expected, this degeneracy breaks the random matrix statistics : this has been shown 
by the computation of the two-point correlation function ([HE])- Moreover it seems 
reasonable to expect that random matrix statistics would be retrieved by glueing some 
star graphs together (just by one edge). Hence star graphs can really be used as a toy 
model for what degeneracy can induce on statistics, and how degeneracy can be broken. 
Their simplicity makes exact results easier to obtain : the trace formula for star graphs 
has been shown to converge under quite reasonable assumptions ([9]). Moreover, star 
graphs may also be considered as a discrete version of Seba billiards (pi)]) : indeed the 
eigenvalues of quantum star graphs and the energy levels of Seba billiards are solutions of 
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similar equations, so that their study can also say something about continuous dynamics, 
not only discrete one. As a step further in the understanding of this model, we will here 
compute the three point correlation function of such graphs. Moreover, it is likely to be a 
useful ingredient for the computation of the two point correlation function for glued star 
graphs. Glueing decreases degeneracy, so that one could expect another intermediate 
statistics (double star graphs have been considered in [II] . but their graphs are not the 
same as here : they have only two vertices). 

Spectral statistics for quantum graphs and the trace formula relating them to 
periodic orbits will be recalled in the first part. The second part will state the 2-point 
correlation function as obtained in [8]. The third part will present the computation of 
the 3-point correlation function. Perspectives regarding glued graphs will be given in 
the conclusion. 

2. Quantum graphs : eigenvalues and trace formula 

We will start by some vocabulary and notations. Let G = (E, V) be a graph with 
a metric structure : to each edge (i,j) G E C V x V is assigned a length i^-, 
such that kj = Ijf, although the graph is supposed to be non oriented, that is 
(i,j) G E =>- (j,i) G E, and ly = we will consider the edges to be oriented : 
(i, j) is different from (j, i), it really describes the edge going from i to j. On each edge 
(i, j), one can thus define a coordinate x such that x = corresponds to the vertex i, 
and x = corresponds to the vertex j. A periodic orbit of length n is a set of n edges 
(pi, . . . ,p n ) such that pi ends where pi + \ starts (as well as p n and pi). A periodic orbit 
is called primitive if it is not the repetition of a shorter periodic orbit. A primitive orbit 
repeated r times is a non-primitive orbit with repetition number r. We will denote by Vj 
the valence of the vertex j, that is the number of its neighbours. On each edge (i, j), one 
is looking for the spectrum of the Laplacian. In other words, one wants to find A and 
tpij such that — ^-pr = \ 2 tpij(x). As one looks for eigenf unctions defined on the whole 
graph, one imposes continuity relations at each vertex, ipijifS) = V'ifc(O)- Moreover, the 
function should have a unique value on a given point, regardless of the sense of the edge 
it belongs to : hence one wants ipij(x) = — x). Finally, one imposes Neumann 

condition on each vertex . ^-|o — 0. It is then a simple exercise to check that the 
eigenvalues A are the solutions of det(J — e~ tXL S) = 0, where S and L are \E\ x \E\ 
matrices : L is diagonal with the length of each edge as diagonal element, and S is 
defined by S^Qj^) = -5 ijk + |-- 

The trace formula as obtained in [1] (a derivation specific to star graphs is given 
in [9]) states that if d(X) = J2 n ^(^ — ^n) i s the spectral density, then d(X) = 
+ i J2 n X] P GP n 7~A p cos(AZp). Here L is the total length of all edges, P n is the 
set of all periodic orbits of period n up to cyclic reordering (that is Po,Pi,P2 and 
Pi,P2,Po are the same orbits), l p is the length of the orbit, r p its repetition number, 
and A p = Yl i=1 S PuPi+1 . 

We will work on star graphs with v + 1 vertices (see figure [T]): these are graphs 
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Figure 1. A star graph 

with V = {0, ..v } and E — {(0, i), (i, 0), 1 < i < v } : t> vertices are all connected to the 
center 0. The ^-matrix elements are 6(0,1), (i,o) — 1 (this corresponds to trivial scattering), 
>%,o),(cm) = — 1 + \ (backscattering) and S^oyoj) = - (normal scattering). The lengths 
will be taken so that they are incommensurate, and that their distribution is peaked 
around 1 : for instance, they can be chosen randomly, uniformly in [1 — l/2t> , 1 + l/2u], 
each length being independent from the other ones. With such a distribution, an orbit 
of period 2k has a length in [2k — k/v,2k + k/v]; such intervals for differents /c's less 
than v do not overlap. The interesting limit will be the limit v tends to infinity : in 
this limit, orbits with the biggest contribution A p will be orbits with a large number of 
backscatterings. 

3. 2-point correlation function 

The 2-point correlation function is defined as R^ix) = (d(\)d(\ — ^f 5 )}. 

The brackets denote a mean value with respect to the A's, that is (/) = 
limA^oo 7^- J A A f(X)d\. Using the trace formula and performing the integral, one gets 

M*) = i + ^E ^AA'Vv cos(^), 
PtP 

where the sum is over the pairs of periodic orbits (p,p'), up to cyclic permutations, of 
lengths l p and l p i and repetition numbers r p and ry. 

A combinatorial analysis of periodic orbits leads ([8]) to the formula 
Rzix) = 1 + /_ °° K{t) exp(— 2inxT)dr, where K is given near r = by 
K{t) = exp(-4T) + J27=2 Em=o fC M T M +i +1 . The C M are defined by 

/ n i -\-k i \ 

fw-f-^V (K+j~l)\(N+j~l)\ r-rj V J with - V j k- 

M V ^) l^k 1 +...+k j +n 1 +...+n j =M (M+j-l)\ 1 U=l (n i +l)!(fe i +l)! ' _ Z^i=l ft » 

and N = Yll=\ n i- 

It is found to be different from the Poisson statistics, since K(r) clearly depends 
on r, and also different from random matrix statistics, since K(0) = 1 and not 0. 
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4. 3-point correlation function 



The 3-point correlation function is defined in a similar way as (^■) 3 (d(X)d(X — -|p)cf(A — 
Using the trace formula and developing, it is equal to Rz(x) + R^y) + #2(2 — 
y) — 2 + Rs(x, y), where 



Rs&y) = (j) 3 Km / A E ^^W^V 
W a-oo2A./_ a < n r v r v ,r p „ 

p,p,p" r r r 



cos(A/ p ) cos((A j—)lp') cos((A -^)/ p »)dA. 

The sum is over triplets of periodic orbits {p,p',p"), up to cyclic permutations, of 
lengths lp,lp',lp" and repetition numbers r p ,r p ',r p ". Performing the integral, one gets 



Rs{x,y) = Y] Iplp'lp" A A A $ cos (^-{yl ,, - xlpi)) 

L p^p„ r p r p' r p" L \ L J 

+^ P -V+V cos (^(v 1 p" ~ xl p')^j + S ip-i P '-i P » cos [^(v 1 p" + xl P')^j 

Let us have a look at the first term enforcing l p n = l p + l p i. The other ones can 
obviously be treated in a similar way. We thus have to deal with 

R\ = ^ £ lMlp + lp>) A p Ap,Ap„ cos ( 2 -^(yl p + (y - s)^ 
^ „ ^ „// r p r p' r p" \ L J 



Now, since the edge lengths are incommensurate, the condition Z p » = l p > + l p implies 
that the orbit p" is formed of the union of the edges of orbit p and orbit p', which we 
will denote by {p") = (p) U (p'). Moreover, since the edge lengths are sharply peaked 
around 1, the length of a periodic orbit of period 2k is nearly 2k and the total length L 
is nearly 2v . 

Thus one gets approximately 



fc./c' ' p period 



2k r p r p' r p" 



p f period 2fc' 
(p")=(p)U(p') 



Following [8J, we will now take as a parameter the number j of distinct edges visited 
by the orbit p" . 



4-1. The j : = 1 case. 

The j = 1 case is a bit special : the orbits p, p' and p" are here all formed of one and 
the same edge, for which there are k) = v choices; the repetition number of such an 
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orbit of period 2k is then r p = k, and the number of backscatterings is also k, thus the 
contribution A p is (—1 + ^) k . Hence, the j = 1 term is 

cy fey \ / o\ 2fc + 2fc' 

^ = ^ E cos ( V (% + k ' iy ~ X)) ) ( _1 + v) 

In the u — > oo limit, putting r = the sum becomes an integral : 

R 1 / -2 [[ drdr'e- 4T e- 4r ' cos(2tt(?/t + (y - x)r')). 
J Jr+ 2 

4-2. The j = 2 case. 

This case is still a bit special, because while p" passes through 2 distinct edges, p and 
p' may still be restricted to one edge only. The orbit p" of period 2k" is here formed 
of 2 different edges, denoted by "a" and "b", for which there are (Z) choices. It can 
be described as a succession of m" packets of a and b. The number of scatterings and 
backscatterings of such an orbit is respectively 2m" (each change of packet contributes) 
and k" — 2m" . Its decomposition in p U p' will obviously also depend on the numbers 
n" of edges a and n' b ' = k" — n" of edges b (we count the edges only when they depart 
from the root). One then needs to know how many different orbits p" there are with 
given k" ,m" ,n",n b , since they will all contribute the same. This number, as explained 
in [8], can be computed as follows : you first divide your a's into m" packets (order 
counts), then your b's into m" packets. Reminding that the number of partitions of an 
integer N into K parts is (^1*), the number of ways to do this is fm"-i)" ^ ne 

orbits should be counted up to cyclic permutation, and the weight of each orbit has to 
be divided by its repetition number r. But one orbit with a given r corresponds to m" jr 
such decompositions. For example, the orbit aabab has m — 2, r — 1, and actually we 
get it twice, since it is the same up to cyclic permutation as abaab. abab has m = 2 
and r = 2, and it is obtained only once. Hence each decomposition gets a 1/m" factor, 
as well as a contribution A" = (-1 + 2)fc"-2m"/2N 2m "_ 

Each such orbit has then to be decomposed into two orbits p and p', composed of 
respectively n a a's and n' a = n" — n a a's (and rib and n' b b's), forming respectively m 
and m' packets. The period of p is 2{n a + nj,). One has to pay attention to the fact that 
n a can be 0, in which case m — 1 but the number of scatterings is then and not 2 
(this does not happen for j > 1 in the case of the 2-point correlation function computed 
in [8] since there are then only two orbits p and p' visiting the same edges). This term 
has to be computed separately. 

Putting all that together, the term R^ 2 where p, p' and p" are all composed of two 
different edges is : 

<-l mm(n a ,n b ) min«,n^) mm«',n' 6 ') 

^' 2 = !(D£ £ £ £ £ £ (< + <)K+<)(n a+ n b ) 

n''n" "o=l "b=i m=l m'=l m"=l 

a ' b / _ 



"'a n,—n. 
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cos(— {(n a + n b )y + (y - x)(n' a + n' b ))) ( -1 + - 

(5:1) (5:Q (:?:!) (-') (TP ^y<»+™'-"> 



One can now perform the v — > oo limit : denoting g* = ^- (the * is either void, ' or "), 
sums over n's turn into integrals over g's , powers of (1 — -) turn into exponentials of 

("-!)! 
(n—m)\ 



g's, terms such as < > n _ turn into (qv ) m 1 . Hence, writing q' = q" — q : 



i£ 2 ~ / dg a dg 6 dg;'dg£(g a + q b M + QbM + 5ft) cos(27r((g a + g fe )y + (y - + $)) 



cxp( 



4(V' + o'MU 3 V l!ME! V K qaqb > V K qaqb > 

m">l v ' m'>\ v 7 m>l v ' 



where the integral is over {q" b > and < g aj ;, < q" b }- Using the modified Bessel 
function h(z) = (z/2) Y^ken k\(k+i)\ ' and denotm g = h{^y/x)/y/x, one gets : 

~ 8 // drdr'e-^e' 4 "' cos(27r(yr + (y - x)t')tt'{t + r') 

J J (M+) 2 

f dg f dg'3((g + q'){r + r' - q - q'))3(q(r - q))3(q'(r' - q')). 
Jo Jo 

Let us now look at the case where one of the orbits (p,p') is composed of only one 
edge; for example, the contribution of the term n a = (and thus n' a = n") is : 

min«X') n'l-l min«,n£) ( <-l \ ( n»-l\ ( <-l\ ( n' b -l\ 

|G)E E E E K + o« + n;)^"- ljU "- in '"" ljl ""- lj 



n",n'/ m"=l m'=l 

a' 6 I II 

n b =n b -n b 



2(m'+m") 



/27T \ / 2\ 2 «'+ n b'- m '- m ") f2\ 

cos ( — Ky + (y - z)« + O) ) ( - 1 + "J ( " J 

Since there are two symbols a and b, and since there are two orbits p and p' that 

— 12 

can be degenerate, the total contribution R 3 ' of the "one degenerate orbit" case is, 
when v tends to infinity : 

R 1 / ~ 8 // drdr / e- 4r e- 4T 'cos(27r(yr + (y - x)t')){t + r') 

J J(R+) 2 

r> [ T dq3(q(r + r' - q))3(q(r' - q)) + r / dgj(g(r + r' - g)p(g(r - g)) 
Jo Jo 



When both orbits p and p' consist of one edge (this problem is specific to the j = 2 
case), the contribution R^ 2 is : 
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R 



1,2 



2 /v 



v 3 \2 



min(n' ' ,n%) 



2ir 



E E «« + <)cos(-«y + <(y-x))) 



<-i \ ( <-i 



<',< m"=l 



"-1 / \ m"-l / 1 



m" n a n b 



9 \ 2n' '+2^'-2m" 
1 + - 



2m" 



drdr'e- 4 ^- 4 "' cos(27r(|/r + (y - x)r')(r + t')3(tt'). 



") 2 



All in all, this gives 



R 1 / = R 1 / + i? 4 ' 2 + j? 4 ' 2 



drdr'e- 4 ^- 4 ^ cos(2vr(yr + (j/ - x)t'){t + r) 



a(rr') + 8rr' f dg f dg'3((g + g') (r + r' - g - q'))3(q(T - g)p(g'(r' - g')) 
Jo Jo 

+8r' / dg3(g(r + r'- q))3(q(r' - q)) + 8r / dq3(q(r + r' - q))3(q(r - q)) 
Jo Jo 



4.3. The j > 2 case. 

The orbit p" is now formed of j edges, denoted by (1, 2, . . . j) for which there are 
choices. We will denote by n" the number of edges i in the orbit p", by m" the number of 
groups of adjacent i, and by n" and m" the corresponding vectors of 1? . For example, 
the orbit p" = 11212332 has j = 3, n" = (3,3,2), and m" = (2,3,1). The period 
of the orbit is 2^y i=1 n" = 2N", the number of scatterings is the number of groups 
M" = YH=i m j an d the number of backscatterings is Yll=i( n j ~ m j) = N" — M" . We 
will denote by Q™" the number of orbits with given j, n" and m", each weighted by 

i/v- 

Such an orbit has to be decomposed into two orbits p and p', consisting respectively 
of rii and v! i = n" — n« edges i, and of and groups of i. Some rii or n[ can be zero, 
but not all of them, and as in the j = 2 case, we will have to consider separately the 
case where all the rii are zero but one. 

To lighten the formulas, we will denote by TJ* T* = TT'T" (where T can be any 
quantity we have defined for p,p' and p"); 1 = f (1, 1, ... , 1); if u and v are vectors, 

U == Yji=i u i-> uV = rii=i w iS M ' = ni=i( u »)' anc ^ u < v means Ui < Vi for all i. 

In the general case where all orbits consist of at least two different edges, we have : 



RY 



n« n -1+ 



k*-M* 



M* 
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{0 < n < n" 
n' = n" — n 
1 < m* < n* 

(Rigorously, we should avoid the case where n or n' is 0, but thanks to the 4r term, its 
contribution will disappear when the sums turn into integrals in the v — > oo limit). 

All we need now is to determine the numbers Q™> ■ The computation is done in [8], 
let's just present the ideas behind. We will count the sequences of (1, 2, . . . j) such that 
there are 7ij i's and m ; groups of i, starting by a group of 1, and not ending by a group 
of 1. Due to cyclic permutations, this is not exactly the same as counting periodic 
orbits, but nearly : for example, the orbit 11212332 corresponds to the n\jr = 2 
sequences 11212332 and 12332112, whereas the orbit 23112311 corresponds to the 
only ni/r = 1 sequence 11231123. Q is then exactly the number of such sequences 
divided by n\. To compute this number, one counts the number of ways to put the i's 
in rrii packets, and then to arrange such packets, starting by 1, keeping the order of the 
groups of a given symbol, not ending with 1, and in such a way that two groups of the 
same symbol are not neighbours. The first step gives a factor JTj=i (m -i)" ^ ne secon d 
step is the most tricky one; it can be evaluated using an exclusion/inclusion principle. 
All in all, this gives 



^n(::0<-^£^(JJn 

i=l x ' Kt<m 1=1 



3 / 1 

rrii - 1 



U - 1 



Now all one has to do is to perform the v — > oo limit. Using Ylm=i 2~2t=i 
Et=i Em =i and Em =i ^Ztyi = xt ~ l ex PX and introducing r = YjUi 3j> one § ets 



*i j =4 e / _rW/ ri d « 

J 1 t)t , it „ A">0 i=1 ^0<q<q», =1 



J P 3 

jTT 1 ' '{t" — t) cos(2it(tx + t"(d — x))) 



n^'--- (r -g^r 1)! ^>— 

Using the identity f >n ^ dq r • -dg ? _i T\U Q?* = / JLi m » ! M+j-i 
developing (q")*' -1 = (q + q')* " -1 , this becomes 

f& ^*JJ +2 drdr'(r + r') cos(2n(yr + (y - x)r')) E E E E 

^' R+ t>l t'>l t">l 0<s<t" 



and 



n 



' ! 1 ' <-i J (T-i)!(r-i)!(T"-r 



(WW*?)' (5 + T — 1)!(T" + T' - S -j - 1)! 

^_2^+T'+T" r 5+T^ r /^T'+T"-S , -i 
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The case where p consists of only one edge gives a factor 



, ) .2N"-M'-M" \ M"+M' 



OS' Qn> [-1 + y 

where J^ 2 denotes a sum over the vectors n", m", m' and over the integer n\ satisfying 

' 1 < m < n'[ 
1 < m" < n" 
1 < m' < n" 
1 < m 1 < n 1 — rti 

Hence, in the limit v — > oo, the contribution of p or p' consisting of one edge only is : 
Rl' j ~ - j — J j ^ drdr'(r + r) cos(27r(yr + (y - x)t')) 

£ (•*;-) (7- - DK7- - di . 2)T , +T „ 

ft i=z (n - 1 - »)i nii«)! nti«?)! < r " + r ' - + s - 

j (r/) T» + T'-t» +s - J+ i r t»-i- s exp( _ 2r ) + T T"+T'-t>{+s- j+ i( T iy>{-i-s exp (-2r / )] 
Finally, the contribution for j > 2 is = + R^ . 

4-4- The 3-point correlation function 

Putting everything back together, the 3-point correlation function can be written 
R 2 (x)+R 2 (y)+R 2 (x-y)-2+ ff R+2 dTdT'(cos(2ir(yT+(y-x)T')+cos(2ir(yT'-x(T+T')) + 
cos(2n(yT + XT'))F(T,T r ), where F = F l + F 2 + F 3 + F A is given by Fi(r,r') = 2e- 4r e - 4r ', 

F 2 (r, r') = e- 4r e- 4T '(r + r') J(rr') + 8r' f dq3{q{r + r' - q))3(q(r' - q)) 

Jo 

+8r rdq3(q(T + r'-q)P(q(r-q)) 
Jo 

+ 8rr' f &q f dq'3((q + q')(r + r' - q - q'))3(q(r - q))3{q\r' - q')) 
Jo Jo 

j ( Si+U-l\ ( t'l-Vt'- Si -l\ 

%,o=e|hoeee e i t « ( J1,£,' ' 

(T - 1)!(T' - 1)!(T" - 1)! / _2) T + T '+ T " r 5+T (r / ) T ' +T "~ s '~-? 



(S* + T — 1)!(T" + T' - S - j - 1)! 
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Figure 2. The function F (line) and its expansion up to order 2 (squares); (b) 
corresponds to r' = 0. 



(T'-1)!(T"-1)! 



-2) 



rpl | rj-iff 



if) 



i\T"+T' -t'l+s-j+l t'l-l- 



(T" + T' — t" + s — j)\ 

+r T" + r-i' 1 ' +s - J+ i (r ^'/-i- exp( _ 2r ') 



exp( 



-2t) 



4-5. Small t,t' expansion. 

The first contribution for the "general" j-th term where no orbit consists of only one edge 
is (tt'^t + r'), and the first contribution when one orbit is degenerate is fYij + r'). 
Keeping only the first terms, and using the expansion of 3(x) = 2 + 4x + 0(x 2 ) when x is 



small, one thus gets F(r, r') = 2 - 6r - 6r' + 16rr' + 8r 2 + 8(r 



A 2 



•or, r 



tt'). Figure 



|2] presents a graph of the function F, together with the first terms of its expansion. 



5. Conclusion and perspectives 

We have obtained here a formula for the 3-point correlation function of star graphs, 
allowing in principle to get its small r, r' expansion at any fixed order s by keeping 
terms up to j = s — 1. One should be able to get the first order of the expansion for 
the n-point correlation function in the same way by computing the j = 1 and j = 2 
contributions. Since the statistics of the eigenvalues are characterized by all the n-point 
correlation functions, this is a small step towards knowing a bit more of this intermediate 
statistics for star graphs. 

This specific result for the 3-point function could also be of some help in computing 
the form factor of two star graphs S\ and S2 glued together by an edge e linking their 
centers. Indeed, the form factor is a sum over pairs of orbits (p, q) visiting the same 
edges. Let us consider all the couples (p, q) corresponding to a given set of edges. To get 
the main term in the small r expansion, it would seem reasonable to group the edges 
corresponding to each graph together, and thus to consider first the couples (p, q) such 
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that p and q can be written p = piEp 2 E and q = q\Eq 2 E, where pi and qi are orbits on 
Si, and E denotes any sequence e...e of edges e. The sum over such p and q should thus 
involve a sum over (pi, q{) and a sum over (p 2 , q 2 ), and the product of the form factors of 
each star graph should appear here. The next term would correspond to decompositions 
p = piEp 2 Ep[Ep 2 and q = q\Eq 2 E. The sum over orbits p, q would then involve a sum 
on each star graph over orbits Pi,p'i,qi such that (<&) = (pi) U (pQ, similar to the sum 
computed in the three-point correlation function. Of course, the exact calculation would 
involve all the n-point correlation functions as well as combinatorial factors to insert the 
£"s, but the first terms of the expansion should already give an insight of what happens. 
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